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Abstract 

Sound propagation in water-saturated granular sediments is known to 
depend on the sediment porosity, but few data in the hterature address 
both the frequency and porosity dependency. To begin to address this 
deficiency, a fiuidized bed technique was used to control the porosity 
of an artificial sediment composed of 265 /xm-diameter glass spheres. 
Time-of-flight measurements and the Fourier phase technique were uti- 
lized to determine the sound speed and attenuation for frequencies from 
300 kHz to 800 kHz, and porosities from 0.37 to 0.43. A Biot-based 
model qualitatively describes the porosity dependence, but a transition 
to a scattering-dominated regime is also observed. 



PACS numbers: 43.30.Ma, 43.35.Bf 
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I. INTRODUCTION 



Numerous models have beenproposed to describe sound propagation in water saturated 
granular materials (see Refs. Q, 0, 4, 5| for a sampling); however, no model yet accurately 
describes both the sound speed and attenuation in these materials across the full range of 
frequencies of interest in underwater acoustics. Measurements of the acoustic properties of 
water saturated sediments have been made under a variety of conditions including in situ 
measurements,-!^ retrieved samples, ^''^ reconstituted samples,-'iS and artificial sediments.— 
There is a large variability in the observed sound speeds and attenuation coefficients of these 
various sediments due to the wide range of porosities, permeabilities, and other physical 
properties they exhibit. Richardson and Jackson have studied the effect of porosity on a 
global scale, compiling numerous measurements obtained across many sediment types.— 
Naturally, the sediment porosity changes are also accompanied by grain size and material 
differences. 

In none of these studies, nor in any other studies known to the authors, was the effect of 
porosity studied in a controlled manner across a range of frequencies, for a single sediment 
sample with uniform grain size distribution. In such an idealized case, one can focus on the 
effect of porosity in isolation, and with sufficient knowledge of most of the sediment physical 
parameters, one can make meaningful comparisons to various models of sound propagation. 
Towards this end, a fluidized bed apparatus was used to prepare specific porosities of a 
spherical glass bead sediment, and ultrasonic time-of-flight measurements were obtained. 
Experimental results were compared to the predictions of the Williams effective density 
fluid model.— The purpose of this brief communication is to describe the porosity-control 
technique as applied to sediment acoustics and to present a set of example measurements. 



II. DESCRIPTION OF EXPERIMENT 



The artificial sediment was composed of 265 j^m ±15 /im diameter soda lime glass 
spheres saturated with degassed distilled water. The beads were sorted using a sieve shaker 
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between 250 /xm and 280 /im fine mesh sieves. The mass of the dry spheres was measured, 
and the spheres were then added to the apparatus shown in Fig. [TJ Degassed distilled water 
was then continuously pumped through the sample for 30 minutes and then allowed to rest 
overnight to dissolve trapped gas bubbles. 

Sediment porosity j3 was controlled by using the fiuidized bed apparatus shown in Fig. [U 
which is similar to that described by Schroter et al.—^ The interior of the sample holder 
was 6 cm x 6 cm square and 30 cm tall, with the sediment height dependent upon porosity, 
but nominally 100 mm. Two parallel polycarbonate walls of the sample holder were 1 mm 
thick and the opposite two were 10 mm thick. Control of the porosity of the monodisperse 
fluid-saturated bead pack is initiated by pumping water at a particular flow rate through 
the distributor (a block of sintered bronze powder, 5 mm thick with approximately 10 /im- 
diameter pores). The beads are carried into suspension, with the overall magnitude of their 
motion governed by the flow rate. The beads achieve steady-state motion and after two 
minutes, the flow is stopped. The beads are then given one minute to settle. The resulting 
porosity is dependent upon the flow rate. An additional procedure was use to achieve 
uniform porosity for (3 < 0.40: A series of flow pulses were introduced, starting at a high 
flow rate but with an incrementally decreasing flow rate. This allowed the bead column to 
settle more uniformly. 

A rectangular excitation pulse with a 10 volt amplitude and a 1.25 /is width was di- 
rected to a pair of broadband ultrasonic transducers (labeled Tx) with center frequencies 
of 500 kHz, which transmitted simultaneously through the water and sediment paths, as 
shown in Fig. [H Another pair of ultrasonic transducers (labeled Rx) with center frequen- 
cies of 550 kHz received the acoustic signals. Commercially-available ultrasonic coupling 
gel was used between the transducer faces and the sample holder walls. A digital storage 
oscilloscope acquired both the source signal and the received signals. In all cases, 512 linear 
time-domain signals were averaged onboard the scope to remove noise. The averaged signals 
were transferred to a computer for analysis. 

To determine sediment porosity, images of the beads were obtained with a digital camera 
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at the same time as the acoustic measurements. A fiducial scale indicating sediment height 
h was placed within the field of view of the camera. The porosity of the sample (3 was 
determined with (3 = 1 — rrib/ {p^hA), where rrib and pb are the dry mass and density of the 
beads, h is the height of the beads within the column, and A is the cross-sectional area of 
the column. A resistance temperature detector, labled RTD in Fig. [H was used to measure 



161] was later 



the temperature during the acoustic measurements. Equation 5.6.8 of Ref. 
used to determine the sound speed cq of the distilled water at the experimental temperature. 

The apparatus and preparation technique described above yielded a uniform porosity 
throughout the sample,— but also resulted in an acoustic path length that was shorter than 
desired. Increasing the cross-sectional area of the apparatus while retaining a homogeneous 
porosity is difficult. Since the goal of this work was to investigate the effect of porosity, we 
chose to achieve a homogeneous sediment at the expense of simplicity in the acoustic data 
analysis. In addition, the porosity control procedure can take as long as several hours for 
the dense sediments. Because the water temperature follows temperature fluctuations in 
the laboratory, there could be signiflcant temperature-depenent sound speed changes in the 
sample before and after the procedure. Hence, the four-transducer apparatus of Fig. [T]was 
used to facilitate contemporaneous post-preparation water-path calibration and sediment- 
path measurements. 

Since the path length was short, on the order of a few wavelengths at the lowest fre- 
quencies, coupling between the sediment sample and the tank walls and transducers was 
found to be signiflcant, and the short path length increased error due to uncertainty in path 
length and path length differences. The coupling effect is demonstrated in Fig. [2l which 
shows the phase of the measured electrical input impedance of one of the receiving acoustic 
transducers when mounted in the apparatus fllled only with water, and when fllled with 
water-saturated beads. There is as much as a thirty degree phase difference between these 
two cases. Further, initial measurements with distilled water revealed that the sound speed 
and attenuation determined by Eqs. [H and [2] applied to the present four-transducer appa- 
ratus included a systematic error due to a small path length difference between the two 



transducer pairs. The effect of these systematic errors was found to be the same order of 
magnitude (a few percent) as the sound speed changes expected from the porosity depen- 
dency. Therefore, the cahbration and error correction procedure described in Section IIIII 
was implemented. Near-field effects were quantified as described in Ref. 17|] and found to 



be insignificant, which was also the case for a similar experiment described in Ref. 



io|. 



III. ANALYSIS 



The Fast Fourier Transforms (FFTs) of the received signals for both the water-paths 
and sediment-paths were calculated. A Fourier phase technique^ was then used to obtain 
a measurement of the apparent sound speed Cgys of the sediment within the experimental 
system. This sound speed was calculated with 
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where cq was the sound speed in the water portion of the fiuidized bed determined by 



Eq. 5.6.8 of Ref. 



16j, I was the length of travel for the acoustic wave within the sample, 
uj = 27r/ was the angular frequency, and A(j) was the difference in phase between the two 
FFT spectra at frequency u. The apparent frequency dependent attenuation of the sediment 



within the system was calculated with 
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where Aq and Agys were the frequency dependent magnitudes of the FFT at frequency u for 
the water and sediment paths, respectively. 

An equivalent circuit transmission line model of the experimental system (Fig. [3^) was 
used to correct for the systematic errors described in Sec. [TTl The model included the 
transducers, coupling gel, tank walls, and the sediment (or water). The passive layers 



and Z. 



in) 



(Fig. [3b) were modeled with 4"^ = p^^^c^^^S^^V^^ = tan [fc(")/(")/2" 
— jZg"''/ sin /c^")/*^") , where /c*^"'^ p*^"), and c^") are the wave number, density, and sound 
speed within the nth component respectively, is the path length in the nth component. 



S**^"^ is the surface area, and is the transformation factor from one physical domain to 
another. The values of the sediment sound speed and attenuation used in the system model 
are denoted by Csed and aged, respectively, and are used to construct both k and c for the 
sediment layer. A Mason equivalent circuit model of a thickness mode piezoelectric resonator 
with a matching layer— was used to model the transducers (Fig. [Sb). The unloaded electrical 
impedance of each transducer was measured and the values for each model element shown 
in Fig. [3]d were found by fitting the model to the impedance measurement, including the use 
of the properties for PZT-IV listed in Table 4.3 of Ref. 201 . Additional material properties 



for the model elements are listed in Table [H Measurements made with distilled degassed 
water occupying both transmission paths allowed for calibration of the transmission path 
lengths in the model. 

To perform the correction, the experimental input signals for the water and sediment 
paths were digitized, FFTs were taken, the resulting spectra were input into the system 
model, and the output spectra were calculated. These modeled output spectra were then 
subjected to the calculations described in Eqs. [T]and[2l which yield modeled system values of 
Csys,m and asys,m- Using Newton-Raphson iteration, Csed and aged were systematically varied 
in the model until the modeled system values Csys,m and asys,m niatched the measured values 
of Csys and agys for each frequency u and for each porosity. The values of Cged and aged that 
achieve this match were taken to be the free-field sound speeds and attenuation coefficients 
of the sediment samples and are reported in Fig. HI 



IV. RESULTS 

In Fig. and b, the frequency dependence of the corrected sound speed and attenu- 
ation measurements are compared to the predictions of the Effective Density Fluid Model 
(EDFM),^^ for the parameters given in Table HTl Negative dispersion is apparent in Fig. 
for all porosities above 550 kHz. Considering just one frequency, say 600 kHz, the varia- 
tion in the measured sound speed as a function of porosity is well-predicted by the EDFM. 
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Measured values of the attenuation in Fig. |Dd are an order of magnitude greater than the 
predictions of the EDFM, however this discrepancy is consistent with other measurements 
reported in the hterature.— liii^i Trendhnes with frequency /^^^ and dependance are also 
shown. At low frequencies, the data are described by an /^/^ dependence as predicted by 
Biot theory. The sharp increase in attenuation to an dependence at frequencies above 
700 kHz is attributable to scattering eSects.—^^^^ 

In Figs. St and d, the porosity dependence of the measured sound speed and attenuation 
are shown and compared to the predictions of the EDFM at three frequencies. Measurements 
at porosities not shown in Fig. and b are shown here. It is apparent that the EDFM does a 
good job of describing the observed relative porosity dependence. Note that all of the sound 
speeds predicted by the EDFM nearly collapse onto a single line. The negative dispersion 
exhibited by the measured sound speed remains apparent. Plotted this way, the sound speed 
monotonically decreases with porosity, while the attenuation is nearly constant. Therefore, 
the porosity appears to be a good indicator of the sound speed but a poor indicator of the 
attenuation in this glass bead sediment. 

V. CONCLUSION 

The porosity- and frequency- dependent sound speed and attenuation in an artificial 
water-saturated glass bead sediment were measured using a fluidized bed apparatus with 
a time-of-flight acoustic measurement and a Fourier phase technique. The relative depen- 
dence of the sound speed and attenuation on porosity is properly predicted by the Biot- 
based EDFM model throughout the experimental frequency range, 300 kHz to 800 kHz, 
and it also quantitatively described the sound speed ratio to within about ±1%. The mea- 
sured frequency dependence reveals negative dispersion at frequencies above 550 kHz. The 
predicted and the measured attenuation is proportional to /^/^ below 700 kHz, but the mea- 
sured attenuation exhibits an proportionality at higher frequencies for all porosities. The 
EDFM underpredicts the absolute value of the attenuation throughout the frequency range, 
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but it is not intended to be accurate at the highest frequencies in this work. The authors 
acknowledge that other models^i^i^ will likely do a better job predicting the dispersion and 
attenuation in this frequency range, but no current model can transition from /^/^ to 
frequency dependence. The present purpose is to report the porosity control technique and 
initial measurements obtained with it, including the observation of the transition into the 
scattering-controlled regime. Measurements with a variety of bead sizes and with natural 
sand grains, and a more thorough model evaluation are currently underway. 
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density (kg/m^) 


sound speed (m/s) 


matching layer 


1200 


2650(1 + 0.01 j) 


coupling gel 


1000 


1500 


sample holder wall 


1200 


1800(1 + 0.01j) 



TABLE I. Parameters used in the transmission line model. The sound speed and density 
of the transducer matching layers and coupling gel were obtained from their respective 
manufacturers. The sample holder wall density was measured, and its sound speed was also 
measured using commercial ultrasonic testing instrumentation. 
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grain density 


2487 kg/m^ 


water density 


998 kg/m^ 


bead bulk mod. 


41.175 GPa 


water sound vel. (m/s) 


1474, 1472, 1471, 1472 


water viscosity 


0.001 Pa-s 


permeability (m^^) 


[5.32, 6.44, 7.50, 8.87] xlO"" 


tortuosity 


1.63, 1.59, 1.57, 1.54 


porosity 


0.376, 0.393, 0.407, 0.423 



TABLE II. Parameters used as inputs to the Williams EDFM model. Grain density was 



16|. 



measured using a Gay-Lussac flask. Water sound speed was given by Eq. 5.6.8 in Ref. 
Tortuosity was calculated using Ref. 2^ with an exponent of 0.5. Permeability was calcu- 
lated using the Kozeny-Carman Equation^^ with a coefficient of 5 and the pore size factor 



given by Ref. 13j. Bead bulk modulus was taken from the manufacturer's specification. 
Porosity was determined by the method in Sec. [Tll When multiple entries appear within 
a cell, they are for the four porosities that are tabulated here, respectively, and that also 
appear in frames (a) and (b). 
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List of Figures 

FIG. 1 A schematic diagram of tiie experimental apparatus is sliown. Solid (dashed) 
lines denote the transmitted (received) electrical signal paths. The glass 
beads are not drawn to scale [15 

FIG. 2 (Color online) The phase of the measured electrical input impedance of one 
of the receiving transducers depicted in Fig. [1] is shown. The blue curve is 
for the transducer mounted in the apparatus when filled with water, and the 
red curve is for water-saturated beads. The phase difference apparent within 
the experimental frequency range (300 kHz through 800 kHz) was found to 
lead to correctable systematic errors in the sound speed measurement. ... [15 

FIG. 3 A schematic of the transmission line model is shown in (a). It is composed of 
a Mason model for the transducers (b) and generic transmission line elements 
(c) for the remaining layers |l6 

FIG. 4 (Color online) Measured sound speed ratios and attenuation coefficients of 
water-saturated glass beads are shown as a function of frequency in (a) and 
(b) and as a function of porosity in (c) and (d). Heavy dashed and solid 
lines correspond to the predictions of the EDFM^ and symbols correspond 
to measurements. Trendlines are shown with light solid black lines in (b) to 
illustrate f^^"^ dependance below 700 kHz and dependance above [l6 
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FIG. 1. A schematic diagram of the experimental apparatus is shown. Sohd (dashed) hnes 
denote the transmitted (received) electrical signal paths. The glass beads are not drawn to 
scale. 




FIG. 2. (Color online) The phase of the measured electrical input impedance of one of 
the receiving transducers depicted in Fig. [1] is shown. The blue curve is for the transducer 
mounted in the apparatus when filled with water, and the red curve is for water-saturated 
beads. The phase difference apparent within the experimental frequency range (300 kHz 
through 800 kHz) was found to lead to correctable systematic errors in the sound speed 
measurement. 
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FIG. 3. A schematic of the transmission hne model is shown in (a). It is composed of 
a Mason model for the transducers (b) and generic transmission line elements (c) for the 
remaining layers. 
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FIG. 4. (Color online) Measured sound speed ratios and attenuation coefficients of water- 
saturated glass beads are shown as a function of frequency in (a) and (b) and as a function 
of porosity in (c) and (d). Heavy dashed and solid lines correspond to the predictions of the 
EDFM— and symbols correspond to measurements. Trendlines are shown with light solid 
black lines in (b) to illustrate /^^^ dependance below 700 kHz and dependance above. 
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